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An  ultraviolet  (UV)  signal  transmission  undergoes  rich  scattering  and  strong  absorption  by  atmospheric  particulates.  We  develop 
a  path  loss  model  for  a  Non-Line-of-Sight  (NLOS)  link.  The  model  is  built  upon  probability  theory  governing  random  migration 
of  photons  in  free  space,  undergoing  scattering,  in  terms  of  angular  direction  and  distance.  The  model  analytically  captures  the 
contributions  of  different  scattering  orders.  Thus  it  relaxes  the  assumptions  of  single  scattering  theory  and  provides  more  realistic 
results.  This  allows  us  to  assess  the  importance  of  high-order  scattering,  such  as  in  a  thick  atmosphere  environment,  where  short 
range  NLOS  UV  communication  is  enhanced  by  hazy  or  foggy  weather.  By  simulation,  it  is  shown  that  the  model  coincides  with  a 
previously  developed  Monte  Carlo  model.  Additional  numerical  examples  are  presented  to  demonstrate  the  effects  of  link  geometry 
and  atmospheric  conditions.  The  results  indicate  the  inherent  tradeoffs  in  beamwidth,  pointing  angles,  range,  absorption,  and 
scattering  and  so  are  valuable  for  NLOS  communication  system  design. 


1.  Introduction 

In  free  space  optical  communication,  the  deep  ultraviolet 
(UV)  spectrum  with  wavelength  200-280  nm  is  regarded  as 
an  appealing  choice  to  overcome  solar  background  radiation 
and  relax  pointing  and  tracking  [3].  High  altitude  ozone 
absorbs  most  solar  radiation  in  this  band,  yielding  negli¬ 
gible  background  noise  at  sea  level  [4].  And,  atmospheric 
scattering  is  very  strong,  enabling  non-line-of-sight  (NLOS) 
communications  where  the  transmitter  is  not  necessarily 
within  the  receiver  field  of  view  (FOV).  However,  in 
addition  to  scattering,  strong  atmospheric  absorption  leads 
to  significant  signal  attenuation  and  so  limits  achievable 
rates. 

These  properties  have  motivated  development  of  UV 
signal  propagation  models  and  communication  systems 
for  military  and  civilian  applications,  for  example,  long- 
range  communication  based  on  high-power  UV  lasers 
since  the  1960s  [5-8].  Recent  progress  in  deep  UV 
light  emitting  diodes  (LEDs)  [9]  and  avalanche  photo¬ 
diodes  (APDs)  [10,  11]  offers  promise  for  deployment 


of  low-cost/power  and  moderate  bandwidth  short-range 
UV  communication  systems  [12,  13],  including  underwa¬ 
ter  communications  [14,  15]  and  sensor  networks  [16, 
17]. 

NLOS  channel  modeling  is  more  complex  than  tra¬ 
ditional  LOS  links.  In  addition  to  wavelength  and  device 
characteristics,  NLOS  path  loss  is  a  function  of  system  geom¬ 
etry,  including  transmitter  (Tx)  beamwidth,  communication 
range,  receiver  (Rx)  FOV,  the  pointing  elevation  angles,  as 
well  as  the  optical  properties  of  the  atmosphere. 

For  simplicity  and  tractable  analysis,  single  scattering 
models  for  NLOS  communication  links  were  developed 
[1,  18],  with  a  corresponding  single  scattering  assumption 
imposed;  that  is,  each  photon  undergoes  only  a  single 
interaction  with  the  atmosphere  before  it  reaches  the 
detector.  Recently  the  model  has  been  further  simplified  for 
improved  analytical  tractability  [19].  However,  the  single 
scattering  assumption  does  not  always  lead  to  accurate  link 
performance  prediction,  especially  as  the  range  increases, 
and  with  beam  pointing  at  low  elevation  angles.  Multiple 
scattering  may  occur  when  the  particle  density  is  large  and/or 
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the  propagation  distance  is  long.  Alternatives  to  the  single 
scattering  model  include  an  empirical  path  loss  model  [20] 
and  a  Monte  Carlo  statistical  path  loss  model  [2].  Those 
models  are  applicable  for  predicting  path  loss  in  a  variety 
of  scenarios  and  are  generally  more  accurate  than  the  single 
scattering  model. 

In  this  paper,  following  the  same  physical  scattering 
law  as  the  Monte  Carlo  statistical  method  [2,  21],  we 
develop  a  stochastic  analytical  NLOS  UV  channel  path  loss 
model.  We  apply  a  stochastic  analytical  technique  [22,  23]  to 
theoretically  derive  the  nth  scattered  signal  energy  collected 
by  the  detector.  The  model  assumes  that  the  photons  are 
stochastically  scattered  and/or  absorbed  by  the  atmospheric 
particles  and  involves  probabilistic  modeling  of  photon 
random  moving  direction,  distance,  energy  loss,  and  receiver 
capture  after  a  specified  number  of  scatterings.  In  order 
to  obtain  the  nth-order  scattered  signal  at  the  receiver,  we 
trace  the  migration  routes  of  a  single  photon  through  the 
medium.  Its  scattering  distance  and  scattering  angles  follow 
certain  probability  distribution  functions  (PDFs)  [24].  The 
propagation  of  a  photon  between  two  consecutive  scatterers 
is  modeled  as  a  single  scattering  event.  The  probability  of 
the  photon  arriving  at  the  receiver  is  a  function  of  the 
n  scattering  events  encountered.  A  similar  technique  has 
been  applied  to  model  multiply  scattered  lidar  returns  in 
cloudy  media  where  significant  scattering  occurs  [22].  To 
account  for  a  divergent  UV  beam  profile,  the  contributions 
of  photons  at  all  possible  directions  within  the  beam  are 
integrated  [25]. 

We  assume  a  homogeneous  atmosphere  with  constant 
scattering  and  absorption  coefficients  and  ignore  atmo¬ 
spheric  turbulence.  This  assumes  that  different  types  of 
particles  are  weU  mixed  and  the  environment  is  stationary. 
Photons  are  scattered  elastically  which  conserves  energy 
but  incurs  energy  loss  during  propagation.  The  detector  is 
assumed  small  in  size  compared  with  the  propagation  space 
and  is  regarded  as  a  point  detector  with  finite  FOV. 

Our  model  offers  an  analytical  formulation  for  NLOS 
scattering  channel  path  loss  and  provides  a  reference  to 
easily  check  other  models  in  a  variety  of  system  parameter 
settings.  We  consider  the  effects  of  different  optical  pointing 
geometries  and  find  that  path  loss  is  relatively  insensitive  to 
the  Tx  beam  angle  but  increases  considerably  with  increased 
Tx/Rx  elevation  angles  and  decreased  Rx  FOV  for  medium 
communication  range.  Not  surprisingly,  numerical  tests 
demonstrate  good  agreement  with  Monte  Carlo  simulation 
[2] .  However,  the  current  analytical  approach  can  distinguish 
the  contributions  of  different  orders  of  multiple  scattering, 
whereas  Monte  Carlo  modeling  can  only  capture  the  total 
scattering  effect.  A  similar  comparison  of  Monte  Carlo 
simulation  and  analytical  methods  has  also  been  described 
by  Lavigne  et  al.  [26] . 

The  model  reveals  when  high-order  scattering  plays 
a  role,  and  thus  offers  more  insight  into  the  channel 
behavior.  Extensive  numerical  tests  show  that  high-order 
scattering  may  play  a  significant  role  for  scenarios  with 
large  baseline  ranges  and  elevation  angles.  Atmospheric 
conditions  also  affect  the  multiple  scattering  contributions. 
A  tenuous  atmosphere  generally  leads  to  higher  path  loss  at 


Figure  1:  NLOS  UV  communication  link  geometry,  depicting 
elevated  transmitter  beam  and  receiver  field-of-view. 


a  small  range  than  a  thick  or  extra  thick  atmosphere  where 
rich  multiple  scattering  enhances  received  signal  strength. 
As  a  by-product,  when  contributions  from  only  single 
scattering  are  retained,  the  corresponding  single  scattering 
path  loss  shows  a  good  match  with  that  predicted  by 
Reilly’s  analytical  single  scattering  model  [  1  ] .  A  disadvantage 
of  the  proposed  model  is  that  it  cannot  provide  pulse 
spreading  information,  while  the  Monte  Carlo  model  can 
[2]. 

The  organization  of  this  paper  is  as  follows.  The 
stochastic  path  loss  model  based  on  a  general  NLOS  UV 
channel  configuration  is  developed  in  Section  2.  It  consists 
of  modeling  photon  direction,  distance  traveled,  probability 
of  arrival  at  the  receiver  after  scattering  n  times,  and  total 
path  loss.  Numerical  analysis  for  the  multiple  scattering 
NLOS  channel  is  carried  out  and  path  loss  performance  is 
compared  with  both  the  single  scattering  model  and  Monte 
Carlo  simulation  in  Section  3.  Channel  characteristics  under 
different  geometric  parameters  are  analyzed  in  Section  4. 
Finally  some  conclusions  are  drawn. 

2.  Stochastic  Path  Loss  Model 

A  NLOS  UV  channel  involves  rich  scattering  and  absorption 
because  of  abundant  suspended  particulates  in  the  atmo¬ 
sphere.  For  NLOS  UV  communication,  scattering  serves  as 
the  vehicle  for  information  exchange  between  the  Tx  and  Rx. 
The  scattered  photons  that  reach  the  receiver  depend  on  the 
link  geometry,  the  atmospheric  optical  properties,  and  their 
random  migration,  as  described  next. 

2.1.  Link  Geometry.  Consider  a  typical  NLOS  communica¬ 
tion  geometry  [19,  20],  as  shown  in  Figure  1.  Denote  the  Tx 
beam  full-width  divergence  by  ai,  the  Rx  FOV  angle  by  a2, 
the  Tx  elevation  angle  by  ^i,  Rx  elevation  angle  by  /J2>  the 
Tx  and  Rx  baseline  separation  by  r,  and  the  distances  of  the 
intersected  (overlap)  volume  V  to  the  Tx  and  Rx  by  ri  and 
r2,  respectively. 

2.2.  Atmospheric  Optical  Properties.  From  the  communica¬ 
tions  viewpoint,  scattering  and  absorption  are  two  dominant 
types  of  photon  interactions  with  the  atmosphere.  To 
describe  a  homogeneous  atmospheric  medium  related  to  UV 
communication,  we  adopt  the  following  coefficients  [1,  26]: 
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the  Rayleigh  (molecular)  scattering  coefficient  Mie 
(aerosol)  scattering  coefficient  absorption  coefficient 
ka,  and  extinction  coefficient  k^.  The  total  scattering  coeffi¬ 
cient  is  defined  as  the  sum  of  the  two  scattering  coefficients 
ks  =  kf^^  +  and  the  extinction  coefficient  is  given 
by  the  sum  of  the  scattering  and  absorption  coefficients  as 

ke  ^5 

The  scattering  phase  function  is  modeled  as  a  combined 
function  of  Rayleigh  and  Mie  scattering  phase  functions 
based  on  the  corresponding  scattering  coefficients  [19,  27]: 


(1) 


scattering,  and  combined  Rayleigh  and  Mie  scattering  has  the 
following  PDFs: 


^  +  (1  -  y)cos20] 


8(1  +2y) 


/““(0)  = 


_(l+^2  -2^cos0) 
0.5(3  cos^e-  1) 

/,  ,  ,n3/2 

(1+g^) 


7 

i^s  rRay 


/<»(0)  =  ^/r(0)+^/rw 


(5) 

(6) 

(7) 


where  ^  =  cos  0  is  defined  from  the  scattering  angle  6.  The 
two-phase  functions  follow  a  generalized  Rayleigh  model  and 
a  generalized  Henyey-Greenstein  function,  respectively: 


p^^y{f,)  = 

l-g2 


3[l  +  3y+(l-y)^^] 
167r(H-2y) 


47r 


1  0.5(3^^-!) 


(2) 


(3) 


where  y,  g,  and  /  are  relevant  model  parameters. 


2.3.  Elementary  Events  for  Photon  Random  Migration.  Gener¬ 
ally,  it  is  impossible  to  predict  with  certainty  the  trajectory  of 
a  photon  that  travels  in  a  medium  of  randomly  distributed 
scattering  and  absorption  particles.  However,  based  on 
single  scattering  theory,  we  can  obtain  the  probabilities  for 
a  photon’s  trajectory  and  calculate  its  multiple  scattering 
arrival  probability  at  the  detector.  The  photon’s  trajectory 
can  be  exactly  described  by  scattering  distance,  scattering 
zenith  angle  9,  and  scattering  azimuth  angle  f.  We  first 
present  the  PDFs  of  these  three  fundamental  elements  for 
our  modeling.  The  distance  between  neighboring  scatters 
can  be  modeled  from  sampling  a  uniform  random  variable 
X  between  zero  and  one  as  r  =  -  In x/ks  [28].  Accordingly 
after  density  transformation  [29],  the  PDF  of  the  scattering 
distance  r  becomes 


where  9  takes  values  in  [0,  tt].  We  consider  a  uniform 
distribution  between  0  and  2n  for  azimuth  scattering 
angle  f  because  there  is  no  azimuth  dependence  in  the 
phase  function  (symmetry  is  assumed).  Thus,  the  azimuth 
scattering  angle  assumes  the  following  PDF: 


2n’ 


(j)  e  [0,27r]. 


(8) 


2.4.  Modeling  of  n  Times  Scattering.  With  those  probabilities, 
we  are  able  to  develop  a  probability  model  for  a  photon 
to  be  received  by  the  detector  after  exactly  n  scatterings. 
Assume  that  a  photon  from  a  UV  source  is  uniformly 
emitted  within  the  beam  divergence  and  centered  at  direction 
angles  do  and  fo,  and  migrates  a  distance  ro  before  the 
first  scattering  happens.  The  solid  angle  within  the  beam 
can  thus  be  modeled  as  a  uniform  random  variable  with 
a  constant  PDF  of  l/flj  where  flj  is  the  source  beam 
solid  angle  as  Oj  =  27r(l  -  cos(ai/2)).  The  probability 
that  it  is  scattered  in  the  infinitesimal  solid  angle  dflo  = 
sin{9o)d9od(l)o  becomes  1/fls sin(0o)d0od0o  [22],  and  the 
probability  that  it  moves  an  incremental  distance  dro  with 
attenuation  is  e“^“''"/r„(ro)dro.  Therefore  the  probability  that 
this  photon  lies  away  from  the  scatter  center  by  ro  and  further 
moves  dro  along  the  infinitesimal  solid  angle  is  the  product 
of  these  two  probabilities: 


g-karo 

dQo  =  —^fro{ro)sin{9o)d9od(l>odro.  (9) 


Mr)  =  k,e-’^‘’-.  (4) 


The  sample  space  for  variable  r  is  [0,  oo  ] .  The  PDFs  for 
zenith  scattering  angle  9  and  azimuth  angle  f  can  be 
further  formulated  based  on  the  scattering  phase  function. 
The  zenith  scattering  angle  under  Rayleigh  scattering,  Mie 


Notice  that  we  adopt  ka  for  exponential  absorption  loss 
within  distance  ro  while  the  scattering  effect  has  been 
absorbed  in  the  scattering  distance  PDF.  After  arriving  at 
the  first  scattering  center,  the  photon  is  scattered  and  then 
continues  to  move  with  attenuation.  The  scattering  center 
is  regarded  as  a  secondary  point  source  emitting  photons 
spatially  following  the  angle  PDFs  described  by  (7)  and  (8). 
From  the  ith  scattering  center  to  the  (i+  l)th  scattering  center 
{i  =  1,2,.. .),  a  similar  small  probability  dQi  conditioned  on 
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all  previous  events  can  be  written  as 

dQi  =  fe,{6i) ((/>,)  sin(0;)d0,d(j();c/ri.  (10) 

The  same  procedure  can  be  successively  applied  for  each  scat¬ 
tering  process.  Assume  that  each  scattering  is  self-governed, 
and  the  distances  and  angles  for  different  scattering  events 
are  conditioned  on  previous  quantities.  Therefore,  the  arrival 
probability  for  a  photon  that  is  scattered  n  times  before 
arriving  at  the  Rx  can  be  derived  based  on  these  transitions. 
Figure  2  shows  the  photon  trajectory  corresponding  to  n 
scatterings. 

After  the  «th  scattering,  we  focus  on  an  infinitesimal 
solid  angle  within  the  receiver  FOV  in  order  for  the  receiver 
to  receive  this  photon.  Define  the  direction  angle  f„  as  the 
angle  between  the  line  connecting  the  receiver  and  the  nth 
scattering  center  and  the  transmitter-receiver  separation  line. 
Since  the  FOV  has  angle  range  of  [(S2  -  0:2/2, /J2  +  0:2/2] 
as  shown  in  Figure  1,  we  can  confine  the  photon  direction 
using  an  indicator  function  I„  which  equals  one  when  the 
condition  (jS2  -  0:2/2)  <  f„  <  {^2  +  0:2/2)  is  satisfied  and  zero 
otherwise.  Therefore,  the  probability  of  the  photon  leaving 
the  Mth  scattering  center  and  reaching  the  detector  becomes 

dQn  =  f6„{^n)UMn)  sin(0„)d0„d(;f>„.  (11) 

We  define  the  photon  arrival  probability  after  n  scatterings 
as  R„;  that  is,  a  photon  undergoes  n  scatterings  and  arrives 
at  the  Rx.  The  photon  trajectory  is  uniquely  described  by 

which  are  scattering 
distance  set,  scattering  zenith  angle  set,  and  scattering 
azimuth  angle  set  for  n  successive  scatterings,  respectively. 
These  variables  are  associated  with  each  other  using  distance 
vectors  rj  for  i  =  0, . . . ,  m,  and  all  involved  distances 
and  scattering  angles  can  be  recursively  obtained  from  the 
previous  distance  vectors  [2,  28].  Finally  we  obtain 


r„  is  needed  because  it  is  a  function  of  the  other  variables. 
Thus,  there  are  a  total  of  3(«  -1 1)  -  1  integration  variables. 

In  our  model,  the  arrival  probability  R„  is  expressed 
through  a  multidimensional  integral,  and  a  closed  form 
is  not  available.  The  dependence  of  r„  and  I„  on  other 
variables  creates  additional  difficulty  in  model  evaluation, 
and  as  n  increases,  the  complexity  increases.  Hence  we  apply 
a  Monte  Carlo  integration  technique,  a  powerful  but  simple 
numerical  integration  method  for  the  approximate  evalua¬ 
tion  of  definite  integrals  [30].  The  idea  is  to  generate  a  large 
number  of  random  sample  points  uniformly  distributed  in 
the  multidimensional  hypercube  D  and  then  calculate  the 
average  value  of  the  integrand  from  the  random  samples 
[22,  30].  In  our  case  we  proceed  as  follows.  There  are  {n  4- 
1)  subspaces  for  3{n  +  1)  -  1  integrals,  represented  by  all 
the  dQi.  We  randomly  generate  sample  points  for  the  first 
n  subspace  integrations,  and  calculate  the  corresponding 
integrands  from  the  generated  sample  points  {ri,6i,<pi)  {i  = 
0,1,..., M  -  1),  denoted  by  vector  rj-  We  then  obtain 
the  location  of  the  nth  scattering  center  and  calculate  the 
quantity  r„  from  all  those  sample  values.  Applying  these 
results,  we  further  calculate  the  index  function  In-  If  it  is 
one,  then  we  randomly  generate  a  sample  point  for  (0„,  0„) 
and  perform  the  integration  over  the  detector  area.  This 
completes  one  sampling  step.  Suppose  that  the  sample  has  a 
sufficiently  large  size  M,  and  denote  the  points  in  the  sample 
by  R 1, ...,  Rm-  Then  the  estimate  for  the  integral  is  given  by 

r  o  ^ 

'’"■]/(«)"  (iw 

where  for  notational  convenience  F{-)  represents  the  inte¬ 
grand  product  of  functions  over  all  dQi,  and  dR  for  all 
differentials  drj-  During  the  process,  we  truncate  the  infinite 
integration  limit  for  all  r,-  (up  to  f  =  «  -  1)  at  a  sufficiently 
large  number  to  ensure  that  the  integrand  has  decayed 
sufficiently  close  to  zero.  We  have  found  generally  that  using 
a  limit  of  several  times  l/fcj  works  well. 

It  is  important  to  clarify  the  difference  between  the 
Monte  Carlo  integration  of  a  definite  integral  as  used  here 
and  the  Monte  Carlo  simulation  of  a  stochastic  process 
[2].  Even  though  both  belong  to  the  Monte  Carlo  family 
and  involve  random  realizations,  Monte  Carlo  integration 
is  a  numerical  method  based  on  the  approximation  of  the 
averaged  deterministic  integrand  function,  whereas  Monte 
Carlo  simulation  is  a  technique  that  relies  on  repeated 
samplings  (numerous  realizations)  of  a  random  process  to 
compute  a  statistical  expectation. 


R«  =  J  J  ■  ■  ■  J  dQo  X  dQi  ■  ■  ■  X  dQn,  (12) 

where  cfQo  is  given  by  (9),  dQi  by  (10),  and  dQn  by  (11). 
The  integration  limits  for  all  variables  cover  their  full  ranges 
except  the  following:  0o  from  (jSi  -  a\/2)  to  (jSi  ■+  ai/2) 
to  ensure  integration  inside  the  source  beam,  and  d„  and 
(/)„  within  the  solid  angle  of  the  receiver  determined  by  the 
receiver  area  and  distance  r„.  Note  that  no  integration  over 


2.5.  Energy  Loss  Model.  To  obtain  the  NFOS  channel  energy 
loss,  we  assume  that  a  UV  source  emits  a  pulse  containing 
N  photons  uniformly  within  the  beam  angle  ai,  and  each 
photon  has  an  energy  e.  So  the  total  transmitted  energy  is 
Et  =  eN.  The  received  energy  after  n  scatterings  can  be 
represented  as 


Erin)  =  EtPn. 


(14) 
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-©-  n  =  1 
-B-  n  =  2 
-0-  n  =  3 

(c)  Received  energy  for  {^1,^2)  =  (60° >  60° ) 


-©-  M  =  1 

-B-  n  =  2 
-0-  n  =  3 

(d)  Received  energy  for  (/?i,  /(2)  =  (90°,  90°) 


Figure  3:  Received  energy  ratio  predicted  by  the  proposed  model. 


The  «th  scattering  path  loss  is  modeled  as 


PL(«) 


Et 

Erin) 


K' 


Note  that  when  n  =  00,  the  energy  represents  the  total 
received  energy  contributed  by  all  scatterings  and  path  loss 
represents  the  actual  link  performance.  If  we  adopt  decibels 
(15)  for  path  loss,  then  it  becomes  10logjo(PL„). 


The  received  total  energy  up  to  n  scatterings  is 


3.  Modeling  Performance 


n  n 

Er,n  =  X^rii)  =  Et^Pi-  (16) 

i=l  i=l 


From  this  result,  the  received  energy  ratio  defined  as  Ey^JEt 
becomes  'Z!i=i  Ph  and  the  corresponding  path  loss  can  be 
expressed  as 


PL  -  A  -  ^ 
”  Er,n  XU  Pi' 


(17) 


In  this  section  we  study  modeling  performance  in  terms 
of  the  received  energy  ratio  and  path  loss  results  for  dif¬ 
ferent  geometry  scenarios.  Through  numerical  simulation, 
multiple  scattering  effects  are  observed  and  analyzed.  These 
results  are  also  compared  with  Reilly’s  analytical  single 
scattering  model  [1]  as  well  as  the  Monte  Carlo  simulated 
multiple  scattering  model  [2],  respectively,  under  the  same 
geometries  and  atmosphere  parameters.  We  demonstrate 
good  agreement  with  [  1  ]  when  single  scattering  is  of  concern, 
and  with  [2]  when  multiple  scatterings  are  considered. 
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(a)  Path  loss  for  (j3i,  P2)  =  (20°,  20°) 


-©-  «  =  1  M  =  3 

-B-  n  =  2  Reilly 

(c)  Path  loss  for  (j3i, /I2)  =  (60°,  60°) 

Figure  4:  Path  losses  predicted  by  the 


-©-  n  =  1  -0-  n  =  3 

-B-  n  =  2  Reilly 

(d)  Path  loss  for  (/?i,  ;S2)  =  (90°,  90°) 
iposed  model  and  Reilly’s  model  [  1  ] . 


We  assume  gas  concentrations  and  optical  features  of 
the  atmosphere  described  by  [7,  Table  II]  for  middle  UV  at 
wavelength  260  nm.  Large  dynamic  ranges  for  the  absorption 
and  Mie  scattering  coefficients  indicate  that  weather  con¬ 
ditions  may  significantly  affect  the  UV  signal  propagation. 
However,  an  explicit  correspondence  between  the  parameter 
settings  and  weather  conditions  is  not  available  in  the 
literature.  Therefore  we  consider  atmosphere  coefficients  for 
typical  tenuous,  thick,  and  extra  thick  atmosphere  conditions 
(corresponding  to  clear,  overcast,  and  foggy),  given  in 
Table  1.  The  tenuous  condition  will  be  adopted  for  all  the 
following  results  unless  stated  otherwise.  This  agrees  with 
the  experimental  measurement  conditions  in  [  1 3  ] .  We  set  the 
geometric  and  model  parameters  as  follows:  (ai,a2)/5iji32) 
=  (17°,  30°,  60°,  60°),  r  ranges  from  10  m  to  1000  m,  y  = 
0.017,  g  =  0.72,  /  =  0.5,  and  the  detector  area  is  1.77  cm^. 
Unless  otherwise  specified,  the  same  parameters  are  used  in 
Section  4. 


To  illustrate  the  importance  of  multiple  scattering, 
consider  the  following  example.  Figure  3  shows  the  ratio 
of  the  accumulated  received  energy  from  n  =  1  (1st 

scattering  order  only)  to  «  =  2  (sum  of  1st  and  2nd- 
order  scatterings),  and  «  =  3  (sum  of  1st,  2nd  and  3rd 
order  scatterings),  all  predicted  via  our  proposed  stochastic 
modeling  for  different  optical  geometries.  We  observe  that, 
for  each  scenario,  the  2nd-order  scattering  contributes  little 
to  the  total  received  energy  for  ranges  within  100  meters, 
but  it  will  increasingly  contribute  more  with  longer  ranges 
and  larger  pointing  angles.  Further,  considering  n  =  3 
does  not  add  much  to  the  received  energy.  For  this  case, 
2nd-order  scattering  is  sufficient  because  3rd  order  scat¬ 
tering  makes  only  a  negligible  contribution  to  the  received 
energy. 

The  received  signal  energy  can  be  transformed  to  the 
link  path  loss  in  decibels,  depicted  for  different  geometries  in 
Figure  4.  We  also  compare  the  proposed  model  with  Reilly’s 


EURASIP  Journal  on  Wireless  Communications  and  Networking 


7 


-©-  Stochastic 
-B-  Simulation 

(c)  Path  loss  for  (/3i, /I2)  =  (60°,  60°) 

Figure  5:  Comparison  of  the  proposed  stochastic 


-©-  Stochastic 
-B-  Simulation 

(d)  Path  loss  for  (j3i, /32)  =(90°,  90°) 
with  the  Monte  Carlo  simulation  model  [2] . 


single  scattering  model.  Both  models  predict  similar  path  loss 
levels  for  M  =  1.  As  expected,  including  multiple  scattering 
reduces  predicted  path  loss.  The  improvement  is  not  obvious 
for  short  range,  consistent  with  Figure  3.  Again,  the  3rd- 
order  scattering  does  not  lead  to  a  significant  improvement 
for  the  received  energy  for  these  cases. 

Due  to  the  fact  that  our  proposed  stochastic  model  has 
the  same  premise  that  underlies  the  Monte  Carlo  simulation, 
we  should  expect  the  similar  path  loss  from  both  models. 
As  an  example.  Figure  5  compares  the  path  loss  results 
generated  from  our  stochastic  model  and  the  Monte  Carlo 
simulation  model.  Four  sets  of  Tx  and  Rx  elevation  angles 
are  included  for  detailed  comparison.  For  each  subplot,  the 
stochastic  model  provides  a  good  match  with  the  simulation 
model  within  the  high  precision  numerical  accuracy  of  the 
simulation.  This  indicates  that  our  proposed  model  is  a 
good  reference  for  evaluating  the  Monte  Carlo  simulation 
model. 


4.  NLOS  UV  Channel  Characteristics 

Since  NLOS  UV  channel  characteristics  are  crucial  to 
communication  system  design,  in  this  section,  we  apply  the 
proposed  multiple  scattering  model  to  further  study  the 
effects  of  different  system  parameters  on  the  channel  path 
loss,  including  link  geometries  and  two  typical  atmosphere 
conditions.  Referring  to  Figure  1,  we  study  angle  sensitivity 
by  varying  one  angle  and  keeping  the  others  fixed. 

Figure  6  illustrates  range-dependent  path  losses  for  (a) 
varying  Tx  elevation  angle  )3i,  (b)  varying  Rx  elevation 
angle  fii,  (c)  varying  Tx  beam  angle  cti,  and  (d)  varying  Rx 
FOV  angle  a2-  Note  that  path  loss  is  not  very  sensitive  to 
the  Tx  beam  angle.  By  contrast,  path  loss  has  the  heaviest 
dependence  on  the  Rx  FOV  angle.  A  wider  FOV  enables 
collection  of  more  scattered  photons  which  results  in  lower 
path  loss.  Path  loss  varies  slightly  with  Rx  elevation  angle 
and  moderately  with  Tx  elevation  angle.  Larger  baseline 
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-©-  20° 
-B-  60° 
^  90° 


(a)  Path  loss  for  different  /Ji 


-©-  17° 
-B-  30° 
^  45° 


-0-  20° 

-B-  60° 

^  90° 

(b)  Path  loss  for  different  p2 


Range  (m) 


-©-  10° 
-a-  30° 
^  45° 


(c)  Path  loss  for  different  ai  (d)  Path  loss  for  different  az 

Figure  6:  Predicted  path  loss  for  different  system  geometries. 


range  tends  to  make  the  angle  effects  more  pronounced.  As 
an  example,  in  the  first  subplot  for  baseline  range  1000  m, 
path  losses  for  three  different  Tx  elevation  angles  show 
differences  up  to  8dB,  while  for  a  shorter  range  of  100  m, 
the  difference  is  within  3dB.  Similar  effects  are  obtained 
for  other  angular  variation.  One  possible  reason  is  that  the 
roles  of  scattering  and  absorption  may  vary  for  different 
geometries.  Larger  propagation  distance  may  bring  about 
more  scattered  photons  to  the  receiver  but  will  also  lead  to 
higher  absorption  loss. 

Up  to  now  we  have  only  considered  tenuous  atmospheric 
conditions.  Next  we  consider  the  impact  of  atmospheric 
conditions  on  path  loss.  Figure  7  depicts  range-dependent 
path  loss  for  tenuous,  thick,  and  extra  thick  atmosphere 


conditions  whose  coefficients  are  given  in  Table  1.  Four  pairs 
of  Tx  and  Rx  elevation  angles  (jSi,fS2)  are  incorporated  for 
performance  analysis.  For  elevation  angle  pair  (20°,  20°), 
path  loss  results  for  extra  thick  atmosphere  are  about  15  dB 
less  than  for  tenuous,  and  7  dB  smaller  than  thick  atmo¬ 
sphere  at  10  m.  As  the  range  increases,  path  loss  increases 
more  dramatically.  The  gaps  between  curves  decrease  until 
a  few  hundred  meters.  Further  range  increase  leads  to  a 
dramatic  path  loss  increase  for  extra  thick  atmosphere. 
As  elevation  angles  increase,  the  lower  path  loss  region 
tends  to  shift  to  a  smaller  range.  For  example,  for  the 
medium  elevation  angle  pairs  (45°,  45°)  and  (60°,  60°),  the 
lower  attenuation  range  shrinks  to  less  than  100  m,  and  for 
elevation  angle  pair  (90°,  90°),  it  reduces  to  less  than  50  m. 
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-0-  Tenuous 
-B-  Thick 

Extra  thick 

(c)  Path  loss  for  (^1,  JS2)  =  (60°,  60°) 


-0-  Tenuous 
-B-  Thick 
-0-  Extra  thick 

(d)  Path  loss  for  (^1, /)2)  =  (90°,  90°) 


Figure  7:  Predicted  path  loss  for  different  atmosphere  conditions. 


Table  1:  Atmosphere  model  parameters. 


Atmosphere 

fc^(km-') 

k“”=(km  ') 

kafkm-i) 

Tenuous 

0.266 

0.284 

0.972 

Thick 

0.292 

1.431 

1.531 

Extra  thick 

1.912 

7.648 

1.684 

In  these  scenarios,  the  scattering  contribution  dominates 
absorption  loss  for  short  ranges.  For  larger  baseline  ranges, 
scattering  loss  increases  as  the  atmosphere  gradually  becomes 
thick  and  thus  causes  large  path  loss.  For  a  given  atmospheric 
condition,  path  loss  increases  by  20  dB  to  40  dB  when 
the  elevation  angle  increases  from  20°  to  90°  at  short 
range.  These  observations  and  predictions  are  helpful  for 
experimental  channel  characterization  and  communication 
system  design. 


In  order  to  further  demonstrate  the  significance  of 
high-order  scatterings,  in  Figure  8  we  plot  the  ratio  of 
the  accumulated  received  energy  from  n  =  I  (Ist-order 
scattering  only)  to  m  =  5  (sum  of  1st  through  5th- 

order  scatterings)  for  an  extra  thick  atmosphere  at  different 
elevation  angle  pairs.  For  each  pointing  scenario,  we  observe 
that  the  Ist-order  scattering  contributes  dominantly  to  the 
total  received  energy  for  short  ranges  up  to  20  m.  As  the 
range  increases,  high-order  scattering  will  contribute  more 
to  the  total  received  energy.  For  example,  the  2nd-,  3rd-, 
and  4th-order  scattering  overwhelm  other  order  scattering 
for  elevation  angle  (90°,  90°)  at  a  range  of  1000  m.  Further, 
we  can  see  that  the  5th-order  scattering  does  not  contribute 
as  much  as  the  lower  order  scatterings.  We  also  depict 
the  corresponding  path  loss  in  Figure  9.  When  considering 
multiple  scattering  contributions,  the  path  loss  decreases  as 
higher  order  scattering  contributions  are  added.  When  com¬ 
pared  with  tenuous  atmosphere,  an  extra  thick  atmosphere 
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Range  (m) 

-©-  n  =  1  n  =  4 

-O-  n  =  2  n  =  5 

^  «  =  3 


Range  (m) 

-©-  «  =  1  n  =  4 

□  n  =  2  n  =  5 

“0“  n  =  3 


(c)  Received  energy  for  {^1,^2)  =  (60°,  60°) 


(d)  Received  energy  for  (^1,^2)  =  (90°,  90°) 


Figure  8:  Received  energy  ratio  for  extra  thick  atmosphere  predicted  by  the  proposed  model. 


yields  significant  high-order  scattering  contributions,  such  as 
3rd-  and  4th-order  scatterings,  which  should  be  accounted 
for  when  predicting  path  loss.  These  results  indicate  the 
remarkable  property  that  for  relatively  shorter  ranges,  NLOS 
UV  communications  links  are  enhanced  in  hazy  or  foggy 
weather. 


5.  Conclusions 

This  paper  proposed  an  analytical  energy  loss  model  for 
NLOS  UV  communication  channels  based  on  the  proba¬ 
bility  theory  of  scattering  and  absorption.  The  model  was 
developed  by  employing  the  PDFs  of  scattering  distance  and 
scattering  angles.  Multiple  scattering  was  incorporated  and 
contributions  of  different  scattering  orders  were  identified. 
The  total  energy  loss  was  modeled  by  summing  over  the 


scattering  order.  The  path  loss  predicted  using  only  the  con¬ 
tribution  from  the  first-order  agreed  with  Reilly’s  analytical 
single  scattering  model,  as  illustrated  with  different  optical 
geometries.  Multiple  scattering  becomes  more  dominant 
for  some  cases,  especially  longer  propagation  distance  and 
larger  elevation  angles.  Our  model  also  provided  results 
that  are  consistent  with  the  Monte  Carlo  simulation  model. 
Channel  characteristics  were  investigated  in  detail,  including 
the  effects  of  varying  system  geometry  and  the  effects  of 
different  atmosphere  conditions. 

Further  study  will  be  conducted  to  develop  an  analytical 
NLOS  UV  channel  impulse  response  model  from  which 
our  path  loss  model  can  be  further  validated  and  channel 
bandwidth  can  be  predicted.  An  analytically  more  tractable 
path  loss  model  is  also  a  future  topic  of  interest,  which 
can  enable  more  intuitive  analysis  for  the  effects  of  system 
geometry  and  parameters. 
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-©-  n  =  l  M  =  4 

-B-  n  =  2  n  =  5 

-0-  ”  =  3 


-©-  «  =  1  n  =  4 

-B-  n  =  2  n  =  5 

n  =  3 


(c)  Path  loss  for  (ji\,^2)  =  (60° ,  60° ) 


(d)  Path  loss  for  (j3i, /I2)  =  (90°,  90°) 


Figure  9:  Predicted  path  loss  for  extra  thick  atmosphere. 


Acknowledgments 

The  authors  gratefully  acknowledge  R.  Drost  for  his  help 
with  the  coding  examples.  This  work  was  supported  in  part 
by  the  U.S.  Army  Research  Office  under  Grant  W911NF- 
09-1-0293  and  the  Army  Research  Laboratory  under  the 
Collaborative  Technology  Alliance  Program,  Cooperative 
Agreement  DAAD 19-0 1-2-00 11. 

References 

[  1  ]  D.  M.  Reilly  and  C.  Warde,  “Temporal  characteristics  of  single¬ 
scatter  radiation,”  Journal  of  the  Optical  Society  of  America  A, 
vol.  69,  no.  3,  pp.  464-470,  1979. 

[2]  H.  Ding,  G.  Chen,  A.  K.  Majumdar,  B.  M.  Sadler,  and 
Z.  Xu,  “Modeling  of  non-line-of-sight  ultraviolet  scattering 
channels  for  communication,”  IEEE  Journal  on  Selected  Areas 
in  Communications,  vol.  27,  no.  9,  pp.  1535-1544,  2009. 


[3]  Z.  Xu  and  B.  M.  Sadler,  “Ultraviolet  communications:  poten¬ 
tial  and  state-of-the-art,”  IEEE  Communications  Magazine,  vol. 
46,  no.  5,  pp.  67-73,  2008. 

[4]  L.  R.  Roller,  Ultraviolet  Radiation,  John  Wiley  &  Sons,  New 
York,  NY,  USA,  2nd  edition,  1965. 

[5]  G.  L.  Harvey,  “A  survey  of  ultraviolet  communication  sys¬ 
tems,”  Tech.  Rep.,  Naval  Research  Laboratory,  Washington, 
DC,  USA,  March  1964. 

[6]  D.  E.  Sunstein,  A  scatter  communication  link  at  ultraviolet 
frequencies,  B.S.  thesis,  MIT,  Cambridge,  Mass,  USA,  1968. 

[7]  D.  M.  Junge,  Non-line-of-sight  electro-optic  laser  communica¬ 
tions  in  the  middle  ultraviolet,  M.S.  thesis.  Naval  Postgraduate 
School,  Monterey,  Calif,  USA,  December  1977. 

[8]  W.  S.  Ross  and  R.  S.  Kennedy,  “An  investigation  of  atmospheric 
optically  scattered  non-line-of-sight  communication  link,” 
type.  Army  Research  Office  Project  Report,  Research  Triangle 
Park,  NC,  USA,  January  1980. 

[9]  M.  Shatalov,  J.  Zhang,  A.  S.  Chitnis,  et  al.,  “Deep  ultraviolet 
light-emitting  diodes  using  quaternary  AllnCaN  multiple 


12 


EURASIP  Journal  on  Wireless  Communications  and  Networking 


quantum  wells,”  IEEE  Journal  on  Selected  Topics  in  Quantum 
Electronics,  vol.  8,  no.  2,  pp.  302-309,  2002. 

[10]  X.  Bai,  D.  Mcintosh,  H.  Liu,  and  J.  C.  Campbell,  “Ultraviolet 
single  photon  detection  with  Geiger-mode  4H-SiC  avalanche 
photodiodes,”  IEEE  Photonics  Technology  Letters,  vol.  19,  no. 
22,  pp.  1822-1824,  2007. 

[11]  S.  C.  Shen,  Y.  Zhang,  D.  Yoo,  et  al.,  “Performance  of  deep 
ultraviolet  GaN  avalanche  photodiodes  grown  by  MOGVD,” 
IEEE  Photonics  Technology  Letters,  vol.  19,  no.  21,  pp.  1744- 
1746,  2007. 

[12]  G.  A.  Shaw,  A.  M.  Siegel,  J.  Model,  and  D.  Greisokh, 
“Recent  progress  in  short-range  ultraviolet  communication,” 
in  Unattended  Ground  Sensor  Technologies  and  Applications 
VII,  Proceedings  of  SPIE,  pp.  214-225,  Orlando,  Fla,  USA, 
March  2005. 

[13]  G.  Chen,  F.  Ahou-Galala,  Z.  Xu,  and  B.  M.  Sadler,  “Experimen¬ 
tal  evaluation  of  LED-based  solar  blind  NLOS  communication 
links,”  Optics  Express,  vol.  16,  no.  19,  pp.  15059-15068,  2008. 

[14]  S.  Arnon  and  D.  Kedar,  “Non-line-of-sight  underwater  optical 
wireless  communication  network,”  Journal  of  the  Optical 
Society  of  America  A,  vol.  26,  no.  3,  pp.  530-539,  2009. 

[15]  D.  Kedard  and  S.  Arnon,  “Subsea  ultraviolet  solar-blind 
broadband  free-space  optics  communication,”  Optical  Engi¬ 
neering,  vol.  48,  no.  4,  Article  ID  046001,  7  pages,  2009. 

[16]  G.  A.  Shaw,  A.  M.  Siegel,  J.  Model,  and  M.  Nischan,  “Field  test¬ 
ing  and  evaluation  of  a  solar-blind  UV  communication  link 
for  unattended  ground  sensors,”  in  Unattended/Unmanned 
Ground,  Ocean,  and  Air  Sensor  Technologies  and  Applications 
VI,  vol.  54 1 7  of  Proceedings  of  SPIE,  pp.  250-26 1 ,  Orlando,  Fla, 
USA,  April  2004. 

[17]  D.  Kedar  and  S.  Arnon,  “Non-line-of-sight  optical  wireless 
sensor  network  operating  in  multiscattering  channel,”  Applied 
Optics,  vol.  45,  no.  33,  pp.  8454-8461,  2006. 

[18]  M.  R.  Luettgen,  J.  FI.  Shapiro,  and  D.  M.  Reilly,  “Non-line-of- 
sight  single-scatter  propagation  model,”  Journal  of  the  Optical 
Society  of  America  A,  vol.  8,  no.  12,  pp.  1964-1972,  1991. 

[19]  Z.  Xu,  FI.  Ding,  B.  M.  Sadler,  and  G.  Chen,  “Analytical 
performance  study  of  solar  blind  non-line-of-sight  ultraviolet 
short-range  communication  links,”  Optics  Letters,  vol.  33,  no. 

16,  pp.  1860-1862,  2008. 

[20]  G.  Ghen,  Z.  Xu,  FI.  Ding,  and  B.  M.  Sadler,  “Path  loss  modeling 
and  performance  trade-off  study  for  short-range  non-line-of- 
sight  ultraviolet  communications,”  Optics  Express,  vol.  17,  no. 
5,  pp.  3929-3940,  2009. 

[21]  A.  N.  Witt,  “Multiple  scattering  in  reflection  nebulae — I:  a 
Monte  Carlo  approach,”  The  Astrophysical  Journal  Supplement 
Series,  vol.  35,  pp.  1-6,  1977. 

[22]  D.  T.  Gillespie,  “Stochastic-analytic  approach  to  the  calcula¬ 
tion  of  multiply  scattered  lidar  returns,”  Journal  of  the  Optical 
Society  of  America  A,  vol.  2,  no.  8,  pp.  1307-1324,  1985. 

[23]  D.  T.  Gillespie,  “Calculation  of  n-scattered  lidar  retrens  for 
large  n  in  an  idealized  cloud,”  Journal  of  the  Optical  Society  of 
America  A,  vol.  4,  no.  3,  pp.  455-464,  1987. 

[24]  A.  FI.  Gandjhakhche,  R.  Nossal,  and  R.  F.  Bonner,  “Scaling 
relationships  for  theories  of  anisotropic  random  walks  applied 
to  tissue  optics,”  Applied  Optics,  vol.  32,  no.  4,  pp.  504-516, 
1993. 

[25]  L.  R.  Bissonnette,  “Multiple-scattering  lidar  equation,”  Applied 
Optics,  vol.  35,  no.  33,  pp.  6449-6465,  1996. 

[26]  C.  Lavigne,  A.  Roblin,  V.  Gutters,  S.  Langlois,  T.  Girasole,  and 
C.  Roze,  “Comparison  of  iterative  and  Monte  Carlo  methods 
for  calculation  of  the  aureole  about  a  point  source  in  the 
Earth’s  atmosphere,”  Applied  Optics,  vol.  38,  no.  30,  pp.  6237- 
6246,  1999. 


[27]  A.  S.  Zachor,  “Aureole  radiance  field  about  a  source  in  a 
scattering-absorbing  medium,”  Applied  Optics,  vol.  17,  no.  12, 
pp.  1911-1922,  1978. 

[28]  S.  A.  Prahl,  M.  Keijzer,  S.  L.  Jacques,  and  A.  J.  Welch,  “A  Monte 
Carlo  model  of  light  propagation  in  tissue,”  in  Dosimetry 
of  Laser  Radiation  in  Medicine  and  Biology,  vol.  IS  5  of 
Proceedings  of  SPIE,  pp.  102-111,  1989. 

[29]  A.  Leon-Garcia,  Probability  and  Random  Processes  for  Electrical 
Engineering,  Addison- Wesley,  Reading,  Mass,  USA,  1994. 

[30]  Wikipedia,  “Monte  Carlo  integration,”  http://en.wikipedia 
.org/wiki/Monte_Carlo -integration. 


Preliminary  call  for  papers 

The  2011  European  Signal  Processing  Conference  (EUSIPCO-2011)  is  the 
nineteenth  in  a  series  of  conferences  promoted  by  the  European  Association  for 
Signal  Processing  (EURASIP,  www.eurasip.orgl.  This  year  edition  will  take  place 
in  Barcelona,  capital  city  of  Catalonia  (Spain),  and  will  be  jointly  organized  by  the 
Centre  Tecnologic  de  Telecomunicacions  de  Catalunya  (CTTC)  and  the 
Universitat  Politecnica  de  Catalunya  (UPC). 

EUSIPCO-2011  will  focus  on  key  aspects  of  signal  processing  theory  and 
applications  as  listed  below.  Acceptance  of  submissions  will  be  based  on  quality, 
relevance  and  originality.  Accepted  papers  will  be  published  in  the  EUSIPCO 
proceedings  and  presented  during  the  conference.  Paper  submissions,  proposals 
for  tutorials  and  proposals  for  special  sessions  are  invited  in,  but  not  limited  to, 
the  following  areas  of  interest. 

Areas  of  Interest 

•  Audio  and  electro-acoustics. 

•  Design,  implementation,  and  applications  of  signal  processing  systems. 

•  Multimedia  signal  processing  and  coding. 

•  Image  and  multidimensional  signal  processing. 

•  Signal  detection  and  estimation. 

•  Sensor  array  and  multi-channel  signal  processing. 

•  Sensor  fusion  in  networked  systems. 

•  Signal  processing  for  communications. 

•  Medical  imaging  and  image  analysis. 

•  Non-stationary,  non-linear  and  non-Gaussian  signal  processing. 

Submissions 

Procedures  to  submit  a  paper  and  proposals  for  special  sessions  and  tutorials  will 
be  detailed  at  www.eusipco2011.org.  Submitted  papers  must  be  camera-ready,  no 
more  than  5  pages  long,  and  conforming  to  the  standard  specified  on  the 
EUSIPCO  2011  web  site.  First  authors  who  are  registered  students  can  participate 
in  the  best  student  paper  competition. 

Important  Deadlines: 


Organizing  Committee 

Honorary  Chair 
Miguel  A.  Lagunas  (CTTC) 

General  Chair 

Ana  1.  Perez-Neira  (UPC) 

General  Vice-Chair 
Carles  Anton-Haro  (CTTC) 

Technical  Program  Chair 
Xavier  Mestre  (CTTC) 

Technical  Program  Co-Chairs 
Javier  Hernando  (UPC) 

Montserrat  Pardas  (UPC) 

Plenary  Talks 
Ferran  Marques  (UPC) 

Yonina  Eldar  (Technion) 

Special  Sessions 

Ignacio  Santaman'a  (Unversidad 

de  Cantabria) 

Mats  Bengtsson  (KTH) 

Finances 

Montserrat  Najar  (UPC) 

Tutorials 
Daniel  P.  Palomar 
(Hong  Kong  UST) 

Beatrice  Pesquet-Popescu  (ENST) 

Publicity 

Stephan  Pfletschinger  (CTTC) 
Monica  Navarro  (CTTC) 

Publications 
Antonio  Pascual  (UPC) 

Carles  Fernandez  (CTTC) 

Industrial  Liaison  &  Exhibits 
Angeliki  Alexiou 
(University  of  Piraeus) 

Albert  Sitja  (CTTC) 

International  Liaison 
Ju  Liu  (Shandong  University-China) 
Jinhong  Yuan  (UNSW-Australia) 
Tamas  Sziranyi  (SZTAKI  -Hungary) 
Rich  Stern  (CMU-USA) 

Ricardo  L.  de  Queiroz  (UNB-Brazil) 


Proposals  for  special  sessions 

Proposals  for  tutorials 

15  Dec  2010 

18  Feb  2011 

Electronic  submission  of  full  papers 

21  Feb  2011 

CTTC^  K; 

Notification  of  acceptance 

23  May  2011 

European  Assoctalion 
for  Signal  Procesmg 

Submission  of  camera-ready  papers 

6  Jun  2011 

de  Catalunya 

Webpage:  www.eusipco2011.org 


